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We have considered the FRW universe in loop quantum cosmology (LQC) model filled 
with the dark matter (perfect fluid with negligible pressure) and the modified Chaplygin 
gas (MCG) type dark energy. We present the Hubble parameter in terms of the observable 
parameters fl m o, &x0 and Hq with the redshift z and the other parameters like A, B, C and 
a. From Stern data set (12 points), we have obtained the bounds of the arbitrary parameters 
by minimizing the \ 2 test. The best-fit values of the parameters are obtained by 66%, 90% 
and 99% confidence levels. Next due to joint analysis with BAO and CMB observations, we 
have also obtained the bounds of the parameters {B, C) by fixing some other parameters a 
and A. From the best fit of distance modulus fi(z) for our theoretical MCG model in LQC, 
we concluded that our model is in agreement with the union2 sample data. 

I. INTRODUCTION 

The combinations of different observations astrophysical data continuously testing the 



theoretical models and the 
arge scale redshift surveys 



sounds of the parameters. Different observations of the SNela 



arge scale redsmit surveys 
13 and WMAP Q 



II, 



5, f|, the measurements of the cosmic microwave background (< 
indicate that our universe is presently expanding with acceleration. 
Standard big bang Cosmology with perfect fluid fails to accommodate the observational fact. 
In Einstein's gravity, the cosmological constant A (which has the equation of state w\ = — 1) 
is a suitable candidate which derive the acceleration, but till now there is no proof of the 
origin of A. Now assume that there is some unknown matter which is responsible for this 



accelerating scenario which has the property that the positive energy density anc 



negative pressure, know as dark energy 
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121 ] . The scalar field or quintessence [13(] is one of 



sufficient 



the most favored candidate of dark energy which produce sufficient negative pressure to drive 
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acceleration in which the potential dominates over the kinetic term. In the present cosmic 
concordance ACDM model the Universe is formed of ~ 26% matter (baryonic + dark matter) 
and ~ 74% of a smooth vacuum energy component. The thermal CMB component contributes 
only about 0.01%, however, its angular power spectrum of temperature anisotropies encode 
important information about the structure formation process and other cosmic observables. 



If we assume a flat universe and further assume that the only energy densities present are 
those corresponding to the non-relativistic dust-like matter and dark energy, then we need to 
know f2 m of the dust-like matter and H{z) to a very high accuracy in order to get a handle on 



fix or wx of the dark energy 



This can be a fairly strong degeneracy for determining 
wx(z) from observations. TONRY data set with the 230 data points (lfj alongwith the 23 



points from Barris et al [17( are valid for z > 0.01. Another data set consists of all the 156 
points in the "gold" sample of Riess et al [4], which includes the latest points observed by HST 
and this covers the redshift range 1 < z < 1.6. In Einstein's gravity and in the flat model of 
the FRW universe, one finds f^A + £}m. — L which are currently favoured strongly by CMBR 



data (for recent WMAP results, see 



10(| ) . In a simple analysis for the most recent RIESS data 



-0.05 
-0.03 



set gives a best-fit value of Q m to be 0.31 ± 0.04. This matches with the value Q, m = 0.29] 
obtained by Riess et al Q|. In comparison, the best-fit £l m for flat models was found to 1 
0.31 ± 0.08 {bj ]. The flat concordance ACDM model remains an excellent fit to the Union2 data 
with the best-fit constant equation of state parameter w = — 0.997^Q054(stat)^Qog2(stat-|-sys 



together) for a flat universe, or w = —1.038 



+0.056/ 



\ +0.093 



359 (stat)_ ( 



(stat+sys together) with curvature 



V 



18j. Cha plyg in gas is the more effective candidate of dark energy with equation of state 
-B/p 19(] with B > 0. It has be en g eneralized to the form p = —B/p a [20j] and thereafter 



2l|]. The MCG best fits with the 3 



modified to the form p = Ap — B / p c 
SDSS data with the choice of parameters A = 0.085 and a = 1.724 
constraints than the previous ones —0.35 < A < 0.025 |23l |. 



year WMAP and the 



221 ] which are improved 



In recent years, loop quantum gravity (LQG) is a outstanding effort to describe the quantum 
effect of our universe. Nowadays several dark energy models are studied in the frame work of 



loop quantum cosmology (LQC). Quintessence and phantom dark energy models |24l. l25j] have 



been studied in the cosmological evolution in LQC. When the Modified Chaplying Gas coupled 
to dark matter in the universe is described in the frame work LQC by Jamil et al j3] who 
resolved the famous cosmic coincidence problem in modern cosmology. In another study [2^] 
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the authors studied the model with an interacting phantom scalar field with an exponential 
potential and deduced that the future singularity appearing in the standard FRW cosmology 
can be avoided by loop quantum effects. Here we assume the FRW universe in LQC model filled 
with the dark matter and the MCG type dark energy. We present the Hubble parameter in 
terms of the observable parameters 0, m , Sl x and Hq with the redshift z. From Stern data set 
(12 points), we obtain the bounds of the arbitrary parameters by minimizing the x 2 test. The 
best-fit values of the parameters are obtained by 66%, 90% and 99% confidence levels. Next due 
to joint analysis with BAO and CMB observations, we also obtain the bounds and the best fit 
values of the parameters (B, C) by fixing some other parameters A and a. From the best fit of 
distance modulus p(z) for our theoretical MCG model in LQC, we concluded that our model is 
in agreement with the union2 sample data. 



II. BASIC EQUATIONS AND SOLUTIONS FOR MCG IN LOOP QUANTUM 

COSMOLOGY 

We consider the flat homogeneous and isotropic universe described by FRW metric, so the 
modified Einstein's field equations in LQC are given by 

and 

* = >+')( 1 -^) < 2 » 

where H is the Hubble parameter defined as H = - with a is the scale factor. Where 
Pc = ie n '^G' 2 n * s ca ^ e< ^ the critical loop quantum density, 7 is the dimensionless Barbero- 
Immirzi parameter. When the energy density of the universe becomes of the same order of 
the critical density p c , this modification becomes dominant and the universe begins to bounce 
and then oscillate forever. Thus the big bang, big rip and other future singularities at semi 
classical regime can be avoided in LQC. Let us note here it has been suggested that 7 ~ 0.2375 
by the black hole thermodynamics in LQC. In LQG, this parameter is fixed by the requirement 
of the validity of Bekenstein-Hawking entropy for the Schwarzschild black hole. The physical 
solutions are allowed only for p < p c . For p = p c , it is called bounce. The maximum value of 
the Hubble factor H is reached for p max = ^ and the maximum value of Hubble factor is • 
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Here p = p m + p x and p = p x , where p m is the density of matter (with vanishing pressure) and 
p x , p x are respectively the energy density and pressure contribution of some dark energy. Now 
we consider the Universe is filled with Modified Chaplygin Gas (MCG) model whose equation 
of state(EOS) is given by 



p x = Ap x 



B 



B > 0,0 < a < 1 



(3) 



We also consider the dark matter and and the dark energy are separately conserved and the 
conservation equations of dark matter and dark energy (MCG) are given by 



Pm + 3Hp m = 



and 



(4) 



(5) 



p x + 3H(p x +p x ) = 
From first conservation equation (4) we have the solution of p m as 

Pm = Pmo(l + z) 3 (6) 
From the conservation equation (5) we have the solution of the energy density as 



Px 



B 



A + 1 



+ C(l + z) 3 ^ 



a+l)(A+l) 



l 

a + 1 



(7) 



where C is the integrating constant, z = - — 1 is the cosmological redshift (choosing ao = 1) 
and the first constant term can be interpreted as the contribution of dark energy. So the above 
equation can be written as 

B 



Px — PxO 



(1 + A)C 



_{1 + A)C + B {l + A)C + B 

x ( 1 + jS )3(o+l)(A+l) 

where p X Q is the present value of the dark energy density. 



(8) 



In the next section, we shall investigate some bounds of the parameters in loop quantum 
cosmology by observational data fitting. The parameters are determined by H{z)-z (Stern), 



BAO and CMB data analysis 



28l432l|. We shall use the x 2 minimization technique (statistical 



data analysis) from Hubble-redshift data set to get the constraints of the parameters of MCG 
model in LQC. 
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III. OBSERVATIONAL DATA ANALYSIS MECHANISM 

From the solution (8) of MCG and defining the dimensionless density parameters O m o = %m 
and fl x o = we have the expression for Hubble parameter H in terms of redshift parameter 
z as follows (8irG = c = 1) 



H(z) = H 



n 



.rO 



B 



(1 + A)C 

+ tz — ^ tt^t — X 



(1 + A)C + B (1 + A)C + B 



. {1 + z) 3(a + l)(A + l)y + ° +nmQ{1 + 



1/2 



Pc 



B 



, (1 + A)C 
(1 + A)C + B + (1 + A)C + B' 



:( i + z) 3(«+i)(A + i)y +a+fimo(1 + 



1/2 



(9) 



From equation (9), we see that the value of H depends on Ho, A, B,C,a, z so the above 
equation can be written as 



H(z) = H E(z) 



(10) 



where 



E{z) 



&x0 



B 



(1 + A)C 



[1 + A)C + B (1 + A)C + B' 



1/2 



B (1 + A)C 

+ A)C + B + (1 + A)C + B X 



1/2 



(11) 



Now E(z) contains four unknown parameters A, B, C and a. Now we will fixing two param- 
eters and by observational data set the relation between the other two parameters will obtain 
and find the bounds of the parameters. 
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z 


H(z) 


a(z) 





73 


± 8 


0.1 


69 


± 12 


0.17 


83 


± 8 


0.27 


77 


± 14 


0.4 


95 


± 17 A 


0.48 


90 


± 60 


0.88 


97 


± 40.4 


0.9 


117 


± 23 


1.3 


168 


± 17.4 


1.43 


177 


± 18.2 


1.53 


140 


± 14 


1.75 


202 


± 40.4 



Table 1: The Hubble parameter H(z) and the standard error cr(z) for different values of redshift 
z. 



A. Analysis with Stern (H(z)-z) Data Set 



Using observed value of Hubble parameter at different redshifts (twelve data points) listed 



in observed Hubble data by 



331 ] we analyze the model. The Hubble parameter H(z) and the 



standard error &(z) for different values of redshift z are given in Table 1. For this purpose we 
first form the x 2 statistics as a sum of standard normal distribution as follows: 



2 _^ (H(z)-H obs (z)) 2 . 

Xstern ~ 2^ ^2U\ ^ ll > 

where H(z) and H b s (z) are theoretical and observational values of Hubble parameter at 
different redshifts respectively and cr(z) is the corresponding error for the particular observation 
given in table 1. Here, H obs is a nuisance parameter and can be safely marginalized. We consider 
the present value of Hubble parameter Hq = 72 ± 8 Kms -1 Mpc" 1 and a fixed prior distribution. 
Here we shall determine the parameters A, B, C and a from minimizing the above distribution 
Xstern- Fbdng t ne two parameters C,a, the relation between the other parameters A,B can 
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be determined by the observational data. The probability distribution function in terms of the 
parameters A,B,C and a can be written as 



L = J e~ Astern p(H )dH (13) 

where P(Hq) is the prior distribution function for Hq. We now plot the graph for different 
confidence levels. In early stage the Chaplygin Gas follow the equation of state P = Ap where 
A < 1. So, as per our theoretical model the two parameters should satisfy the two inequalities 
A < 1 and B > 0. Now our best fit analysis with Stern observational data support the theoretical 
range of the parameters. The 66% (solid, blue), 90% (dashed, red) and 99% (dashed, black) 
contours are plotted in figures 1, 2 and 4 for a = 0.5 and A = 1, 1/3, —1/3. The best fit values 
of B and C are tabulated in Table 2. 



A 


B 


C 




1 


0.904 


0.565 


10.828 


i 

3 


0.561 


0.778 


8.230 


1 

3 


0.849 


0.599 


7.057 



Table 2: H{z)-z (Stern): The best fit values of B, C and the minimum values of x 2 for different 
values of A. 



B. Joint Analysis with Stern + BAO Data Sets 



The method of joint 
has been proposed by 



analysis, the Baryon Acoustic Oscillation (BAO) peak parameter value 



34| and we shall use their approach. Sloan Digital Sky Survey (SDSS) 
survey is one of the first redshift survey by which the BAO signal has been directly detected at 
a scale ~ 100 MPc. The said analysis is actually the combination of angular diameter distance 
and Hubble parameter at that redshift. This analysis is independent of the measurement of 
Hq and not containing any particular dark energy. Here we examine the parameters B and 
C for Chaplygin gas model from the measurements of the BAO peak for low redshift (with 
range < z < 0.35) using standard x 2 analysis. The error is corresponding to the standard 
deviation, where we consider Gaussian distribution. Low-redshift distance measurements is a 
lightly dependent on different cosmological parameters, the equation of state of dark energy and 



s 



H(z)-z data (Stern) A=l 



I < / 



H(z)-zdata (Stern) A=l/3 



\ w 



H(z)-z data (Stern) A=-l/3 



2.0 0.0 



2.0 



Fig.l 



Fig.2 



Fig.3 



Figs. 1-3 show that the variation of C with B for a = 0.5, fiom = 0.29, Q^o = 0.72 with 
A = 1, 1/3, —1/3 respectively for different confidence levels. The 66% (solid, blue), 90% (dashed, red) 
and 99% (dashed, black) contours are plotted in these figures for the H(z)-z (Stern) analysis. 



have the ability to measure the Hubble constant Hq directly. The BAO peak parameter may be 
denned by 

A= ^5 (1 fjM 2/3 (14) 

Here E(z) = H(z)/Hq is the normalized Hubble parameter, the redshift z\ = 0.35 is the typical 
redshift of the SDSS sample and the integration term is the dimensionless comoving distance to 
the to the redshift z\ The value of the parameter A for the flat model of the universe is given by 



A = 0.469 ±0.017 using SDSS data 



34 ] from luminous red galaxies survey. Now the x 2 function 



for the BAO measurement can be written as 



2 _ (A - 0.469) 2 
Xbao ~ (0.017) 2 (15) 



Now the total joint data analysis (Stern+BAO) for the x function may be defined by 



Xtotal — Xstern + XBAO (16) 

According to our analysis the joint scheme gives the best fit values of B and C in Table 3. 
Finally we draw the contours A vs B for the 66% (solid, blue), 90% (dashed, red) and 99% 
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H(z)-z data (Stern)+BAO A= 



0-0 |- i , ,• , I' 

0.0 0.2 0.4 0.6 0.8 1.0 



Fig.4 
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Fig.5 



H(z)-z data (Stern)+BAO A=-l/3 



Fig.6 



The contours are drawn for 66% (solid, blue), 90% (dashed, red) and 99% (dashed, black) confidence 
levels in figs. 4-6 which show the variations of C against B for a = 0.5, fiom = 0.29, ft x o = 0.72 with 
A = 1 , i/3, —1/3 respectively for the H(z)-z+BAO joint analysis. 



(dashed, black) confidence limits depicted in figures 4 — 6 for a = 0.5 and A = 1, 1/3, —1/3. 



A 


B 


C 


y 2 ■ 


1 


0.735 


0.610 


827.909 


l 

3 


0.921 


0.735 


768.499 


1 

3 


0.585 


0.998 


767.440 



Table 3: H(z)-z (Stern) + BAO : The best fit values of B, C and the minimum values of x 2 
for different values of A. 

C. Joint Analysis with Stern + BAO + CMB Data Sets 



One interesting geometrical probe of dark energy can be determined by the angular scale of 
the first acoustic peak through angular scale of the sound horizon at the surface of last scattering 
which is encoded in the CMB power spectrum Cosmic Microwave Background (CMB) shift 
parameter is denned by [351-4371] . It is not sensitive with respect to perturbations but are suitable 
to constrain model parameter. The CMB power spectrum first peak is the shift parameter which 
is given by 
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H(z)-z data (Stern)+BAO+CMB A= 



0.0 0.5 



H{z)-z data (Stern)+BAO+CMB A=l/3 



H(z)-z data (Stern)+BAO+CMB A=-l/3 



U 1.0 - I i 



0.0 |- 

0.0 0.2 0.4 0.6 0.8 1.0 



4 5 6 



Fig.7 



Fig. 



Fig.9 



The contours are drawn for 66% (solid, blue), 90% (dashed, red) and 99% (dashed, black) confidence 
levels in figs. 7-9 which show the variations of C against B for 
a = 0.5, C = 0.01, D, 0m = 0.29, n xQ = 0.72 with A = 1, 1/3, -1/3 respectively for the 

H (z)-2+BA0+CMB analysis. 



1Z 



22 dz 



(17) 



where Z2 is the value of redshift at the last scattering surface. From WMAP7 data of the 
work of Komatsu et al 38| the value of the parameter has obtained as TZ = 1.726 ± 0.018 at the 
redshift z = 1091.3. Now the x 2 function for the CMB measurement can be written as 



2 _ {TZ - 1.726)" 

XCMB ~ 



(18) 



(0.018) 2 

Now when we consider three cosmological tests together, the total joint data analysis 
(Stern+BAO+CMB) for the \ 2 function may be defined by 



XTOTAL — X Stern + XbAO + XCMB (19) 

Now the best fit values of B and C for joint analysis of BAO and CMB with Stern observational 
data support the theoretical range of the parameters given in Table 4. The 66% (solid, blue), 
90% (dashed, red) and 99% (dashed, black) contours are plotted in figures 7-9 for a = 0.5 and 
A = 1,1/3,-1/3. 
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A 


B 


C 


x 2 ■ 


1 


0.735 


0.457 


10022.6 


1 

3 


0.692 


0.694 


9963.5 


1 

3 


0.731 


0.599 


9962.11 



Table 4: H{z)-z (Stern) + BAO + CMB : The best fit values of B, C and the minimum values 
of x 2 f° r different values of A. 



D. Redshift-Magnitude Observations from Supernovae Type la 



The Supernova Type la experiments provided the main evidence for the existence of dark 
energy. Since 1995, two teams of High-z Supernova Search and the Supernova Cosmology 
Project have discovered several type la supernovas at the high redshifts [l - 4 1 . The observations 



directly measure the distance modulus of a Supernovae and its redshift z 



39 



40J]. Now, take 



recent observationa 
the Union2 sample 



data, including SNe la which consists of 557 data points and belongs to 



From the observations, the luminosity distance cIl(z) determines the dark energy density and 
is defined by 



d L {z) = (l + z)H £^L (20) 

and the distance modulus (distance between absolute and apparent luminosity of a distance 
object) for Supernovas is given by 



(i(z) = 5 log 



10 



d L (z)/H 



+ 25 



(21) 



1 MPc 

The best fit of distance modulus as a function fj,(z) of redshift z for our theoretical model 
and the Supernova Type la Union2 sample are drawn in figure 10 for our best fit values of a, 
A, B and C. From the curves, we see that the theoretical MCG model in LQC is in agreement 
with the union2 sample data. 
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Fig. 10 



In fig. 10, u(z) vs z is plotted for our model (solid line) and the Union2 sample (dotted points). 



IV. DISCUSSIONS 



Modified Chaplygin gas (MCG) is one of the candidate of unified dark matter-dark energy 
model. We have considered the FRW universe in loop quantum cosmology (LQC) model 
filled with the dark matter (perfect fluid with negligible pressure) and the modified Chaplygin 
gas (MCG) type dark energy. We present the Hubble parameter in terms of the observable 
parameters fi m o, an d Hq with the redshift z and the other parameters like A, B, C and a. 
We have chosen the observed values of Q m o = 0.28, £I x q = 0.72 and Hq = 72 Kms -1 Mpc -1 . 
From Stern data set (12 points), we have obtained the bounds of the arbitrary parameters by 
minimizing the x 2 test. Next due to joint analysis of BAO and CMB observations, we have 
also obtained the best fit values and the bounds of the parameters (B, C) by fixing some other 
parameters A = 1, 1/3, — 1/3 and a = 0.5. The best-fit values and bounds of the parameters are 
obtained by 66%, 90% and 99% confidence levels are shown in figures 1-9 for Stern, Stern+BAO 
and Stern+BAO+CMB analysis. The distance modulus fi(z) against redshift z has been 
drawn in figure 10 for our theoretical model of the MCG in LQC for the best fit values of the 
parameters and the observed SNe la Union2 data sample. So our predicted theoretical MCG 
model in LQC permitted the observational data sets. The observations do in fact severely 
constrain the nature of allowed composition of matter-energy by constraining the range of the 
values of the parameters for a physically viable MCG in LQC model. We have checked that 
when p c is large, the best fit values of the parameters and other results of LQC model in MCG 
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coincide with the results of the ref. [29] in Einstein's gravity. When p c is small, the best fit 
values of the parameters and the bounds of parameters spaces in different confidence levels in 
LQC distinguished from Einstein's gravity for MCG dark energy model. Also, in particular, if 
we consider the generalized Chaplygin gas (A = 0), the best fit value of critical Barbero-Immirzi 
parameter 7 is 0.2486, where we have assumed the values of other parameters a = 0.5 and 
B = 0.561 for our convenience. In summary, the conclusion of this discussion suggests that 
even though the effect that quantum aspect of gravity have on the CMB are small, cosmological 
observation can put upper bounds on the magnitude of the correction coming from quantum 
gravity that may be closer to the theoretical expectation than what one would expect. 
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